System-reservoir coupling derived by Maxwell's boundary conditions from weak to 

ultrastrong light-matter coupling regime 
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In the standard theory of cavity quantum electrodynamics (QED), coupling between photons 
inside and outside a cavity (cavity system and photonic reservoir) is given conserving the total 
number of photons. However, when the cavity photons (ultra) strongly interact with atoms or 
excitations in matters, the system-reservoir coupling must be determined from a more fundamental 
viewpoint. Based on the Maxwell's boundary conditions in the QED theory for dielectric media, 
we derive the quantum Langevin equation and input-output relation, in which the total number of 
polaritons (not photons) inside the cavity and photons outside is conserved. 

PACS numbers: 03.65. Yz,42.50.Pq,71.36.+c 



The dissipation has long been discussed as an in- 
evitable phenomenon in most systems, and the light is a 
typical target in the study of such open systems. When 
the light is confined in a cavity consisting of mirrors or 
distributed Bragg reflectors, discrete cavity modes are 
well identified, while they have finite broadenings due to 
the loss through the mirrors. In the standard theory of 
quantum optics [l|, Q , coupling between the cavity modes 
and external photonic field is usually supposed as 



h: 



cavity 
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Here, a m is the annihilation operator of a photon in ru- 
th cavity mode, d(w) is the one outside the cavity with 
frequency uj, and K m (uj) is the dissipation rate of the ru- 
th, mode. This expression has successfully reproduced 
experimental results, even when the cavity photons in- 
teract with atoms or excitations in matters. However, 
in the ultrastrong light-matter coupling regime, where 
the Rabi splitting frequency g is comparable to or larger 
than the transition frequency uj cx of excitation in matter 
we encounter a problem of the treatment of the 
system-reservoir coupling [TBI, EH- This is because the ro- 
tating wave approximation (RWA) cannot be applied on 
the light-matter coupling, and the total number of pho- 
tons and excitations is no longer conserved. Then, while 
the number of photons inside and outside the cavity is 
conserved in Eq. (JT|) , we should reconsider the validity of 
this expression carefully. 

In the ultrastrong light-matter coupling regime, even in 
the ground state of the coupled (polariton) system, there 
are virtual photons represented as a "squeezed" vacuum 
state [Ij]. As pointed out by Glauber and Lewenstein 
[ri| , such virtual photons exist even in a simple dielectric 
medium, and its "squeezing" character is different from 
that of the squeezed light in vacuum. Whereas the elec- 
tromagnetic fields are certainly sub- and super-fluctuant 
in dielectrics, such a "squeezed" quantum fluctuation re- 
covers to the one of the coherent or vacuum state when 



the fields escape from the dielectrics to the vacuum. In 
this way, even in the ultrastrong light-matter coupling 
regime, the polaritons simply represent the electromag- 
netic fields in dielectrics, and we cannot generate non- 
classical light outside the cavity at least in the linear 
optical process with classical inputs. Certainly, classi- 
cal outputs are obtained by classical inputs at least in 
the approach of Langevin equations 0, G3] ■ However, 
when we simply suppose Eq. ([T]), we encounter a deli- 
cate problem: Even if the outside is in the vacuum, since 
the virtual photons inside the cavity can escape to the 
outside, the polariton system is inevitably excited flij ]. 

Instead of supposing Eq. ([I]) , the excitation of polari- 
tons by the vacuum can be avoided by considering 
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a\u)pj -p]a{u) , (2) 



where pj annihilates a polariton in j-th eigen-mode and 
Kj(oj) is its dissipation rate. Although this expression 
can be obtaine d by applying the RWA on Eq. ([1]) in the 
polariton base [15|, we will show that expression (j2J is 
rather natural based on the quantum electrodynamics 
(QED) theory for dielectric media fl8l - [2~4j |. then the po- 
lariton system is in principle not excited by the vacuum 
even in the absence of the RWA. We will also compare 
the dissipation rate derived from Eq. ([T]) with the RWA 
and that from Eq. ^ , which show qualitatively different 
behaviors with the increase of light-matter coupling into 
the ultrastrong regime. 

First, we consider a loss-less homogeneous dielectric 
medium, in which photons interact with infinite-mass ex- 
citations 0, [25[ . The Hamiltonian is represented as 

00 

H p = ^2 \hc\k\a) k a k + huj c J)\b k 
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Here, d k and b k are annihilation operators of a photon 
and an excitation with wavenumber k in z direction, re- 
spectively, and satisfy = = 5k, k'- Using 
the photon operator, the vector potential is represented 

MS 



Tl((£>) 



A(z) = 



k—~oo 



2e c\k\SL 



a k 



(4) 



where c is the speed of light in vacuum, £o is the vacuum 
permittivity, S is the area in x — y plane, and L is the 
length in z direction. w ox is the frequency of excitations, 
gk is the light-matter coupling strength, and the coeffi- 
cient of the last term is Dk = g/ £ 2 /w cx . Introducing the 
annihilation operators of lower and upper ( j = L and U) 
polaritons as 

Pj.k = Wjk&k + Xjkh + Ujk^_ k + z jk tf_ k , (5) 

we can diagonalize Eq. © 0, [25[ : 



*p= E E 

j=L,U k— — oo 



(6) 



The eigen-frequencies Uj k and coefficients 
{wjk,Xj k ,yjk,Zjk} are determined by [p jk ,H p ] = 



hujjkpjk and \pjk,Pfk'} 
c 2 fc 2 



Sjj'Sk, k ' 



and w = ujjk satisfies 



47r/3w e 



- (w + i0+) 2 ' 



(7) 



Here, e p (uA is the dielectric function of the polariton 
medium 2jJ, |25| , and we suppose the coefficient Airfi = 
4:c\k\g k 2 /uj cx 3 does not depend on k for simplicity. By 
using the polariton operator, the positive-frequency com- 
ponent of vector potential (jlj is rewritten as |26j 



^) = E E 

j=L,U k— — oc 



hv g (uj jk ) ikz 



2s cujj k n p (ujj k )SL 



Here, n p (w) = y/e p (uj) is the refractive index, and 

, , duo c 



dk n p (uj) + ui{dn p /doj) 



(8) 



(9) 



is the group velocity. In this way, polaritons represent 
the eigen-states in dielectrics even in the ultrastrong 
light-matter coupling regime in 6 > 1. As discussed in 
Ref . [Tj} , the quantum fluctuation of the electromagnetic 
fields is modulated in dielectrics, which can also be ob- 
tained from Eq. © 

Next, in order to introduce boundaries of the polariton 
system, we employ the QED theory for inhomogeneous 
media [2CH24| . We simply consider an one-dimensional 
(ID) system with dielectric function e(z,uj) depending 
on position z and frequency u satisfying the Kramers- 
Kronig relation, and the electric and magnetic fields are 
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FIG. 1. Sketch of the considered one-dimensional system with 
dielectric function in Eq. (|13l) . 

in the x — y plane. The positive-frequency component 
A + (z,u) of vector potential in this system obeys 



dz 2 



-£(z,u>) 



A+(z,w) = -hoJ n {z,lj). (10) 



This has exactly the same form as the wave equation 
derived from the Maxwell equations. Here, fio is the vac- 
uum permeability, and the quantum fluctuation of the 
electromagnetic fields is described by the noise current 
operator Jn(z,w) satisfying 



[Jn(^,w), Jn(z', w')^] = 5(l}— l}')5(z— z 1 ) Im[e(z, a;)]. 

(11) 

The electric and magnetic fields are represented as 
E(z, t) = -(d/dt)A(z, t) and B(z, t) = (d/dz)A(z, t), re- 
spectively, and their positive-frequency components sat- 
isfy the Maxwell equations 



1 d - 
Mo oz 



iljB + (z, lu), 
—iujD + (z, to). 



(12a) 
(12b) 



The latter equation is equivalent with Eq. (|10[) . and 
the displacement field includes the noise current as 
D + (z,uj) = eoe(z,uj) E + ( z,uj) + (i/w) Jn(z, w). Based 
on this formalism [20h24|. the positive- and negative- 
frequency components never mix with each other at least 
in the linear optical process. Then, for polariton system 
confined in an optical cavity, the annihilation operator 
■pj of polariton couples with a(ui) of photon outside the 
cavity, and they never couple with creation operators [pt 
and cr (w)] in the linear optical process. Next, we explic- 
itly consider a cavity system, and derive the quantum 
Langevin equation and input-output relation from the 
above Maxwell equations. 
As discussed in Ref. 



27] . we consider a cavity system 



shown in Fig. [T] where the dielectric function is given as 



e(z, uj) = r]{ui)5{z) + 



1 z < 

£p(oj) < Z < I 



(13) 
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There is a perfect mirror at z — I, and the other mirror 
is placed at z = 0. T)(w) determines the transparency 
between the cavity and the outside. For the solution 
to Eq. (|10p . we suppose the vector potential in the two 
regions as depicted in Fig. [TJ 

,{to)e l{u/c)z + i+ < He~ i ^/ c ) z , 



AJ (z,u) = A 
A+{z,uj) = A+(uj)sm[k p {uj)(l - z)], 



(14a) 
(14b) 

where k p {ui) — n p (uj)uj / c. Since the electric field is com- 
pletely zero at the boundary with the perfect mirror, the 
intra-cavity mode has no amplitude at z = I. Here, as 
discussed in Ref. [20l - |24j ] . the incoming field Ai> can be 
simply derived as 



(15) 



where a> (w) is defined as 
a>(w) = 



dz e~ l{ul/c)z Mz,uj). (16) 



From Eq. (jlip . this operator satisfies [a>(w), at,(a;')] = 
8(u) — u}'), and it corresponds to the annihilation operator 
of an incoming photon. 

Next, we consider boundary conditions determining 
Aj^ui) and AJiuj). From the continuous condition at 
z = derived from Eq. (|12a[) . we get 

Af^cj) + Af < (u) = A+(uj)sin[k p (uj)l}. (17a) 

From the integral form of Eq. (|12b[) , we also get 

[j4i>(w) - Af K (uj)] - in p (u>)A% (ui) cos[fc p (w)Z] 

= -iA{uj)A+(u)sm[k p (uj)ll (17b) 

where A(ui) — rj(uj)uj/c. Solving these Maxwell's bound- 
ary conditions (MBCs), we get 



2A+ > {u 



[1 — ivl(o;)] sin[fc p (w)Z] + in p (w) cos[fc p (w)Z] ' 

(18) 

and Af^ui) is expressed by Eq. (|17a[) . They are also de- 
rived from Eq. (ITU)) by using the Green's function [26j ]. 
As seen in Eq. (|18p . we can find that resonances are ob- 
tained at lo = Q m satisfying 

t&n[n p ((2 m )(2 m l/c] = n p {£2 m ) / A(Q m ) (19) 

and the frequency broadening is proportional to 
A(f2 m )~ 2 . Here, as discussed in Ref. [27J, we consider 
the good cavity limit A{Q m ) 3> n p (f2 m ), and we suppose 
that A(lo) varies only slightly around lj — fl m for sim- 
plicity. Under these assumptions, around the resonance 
uj oj J? m , Eqs. (IT5|) and (I17ap are approximately rewritten 
as 



A+{u) 



<2v g (Q m ) \y/ KMBc(^m) 



(20a) 



(20b) 

where the dissipation rate kmbc(^) is defined as 

KmbcH " AH 5 * ' ( } 

In the semi-infinite region 1, A^^ui) is expressed by the 
annihilation operator ai<(w) of an outgoing photon as 



A+ < (u>) = 



iTTEnCwS 



oi<(w). 



(22) 



On the other hand, in the finite region 2, in the simi- 
lar way for deriving Eq. ((5J, the expression of A^{z) is 
represented in the good cavity limit: 



j—L,U m 



(23) 

where fc m = m7r/Z, i? im = uij,k=k m and [pj m >p]/ m '] = 
5jj'8 m , m >. Then, from Eqs. fTB]), ([2 % and 1(25 )1 . Eqs. flU 
are rewritten in the good cavity limit as 



Pjm(u) = ~ 7j -— -— — — a in (w), (24a) 



j=L,U m=l 

, (24b) 

where the input and output operators are defined as 
Sin(w) = ai > (w)/v / 27r and a out (w) = ai < (aj)/v / 27r. 



Eqs. (|24| have certainly the same form as the quan- 
tum Langevin equation and input-output relation derived 
from Eq. ([2]) by the well-known treatment in quantum op- 
tics [l|, (2| . In this way, we should suppose Eq. ([2]) for the 
coupling between the polariton system and the photonic 
reservoir, instead of Eq. ([TJ with the RWA. 

To check the consistency with the well-known discus- 
sions [H, 0], we suppose the simplified case where the 
resonance frequency of the lowest bare cavity mode is 
tuned to w cx , and we calculate the dissipation rates de- 
rived from Eq. (|2"T)) . For simplicity, we suppose A(w) is 
constant in the frequency range of interest. In the good 
cavity limit A 3> 1, the cavity length is determined as 
Iq ~ ttc/uj cx satisfying tan(w ex ^/c) -C 1. Then, from 
Eq. (f2~T)l . the dissipation rate of bare cavity mode is ob- 
tained as n — 2c/A 2 l . From Eq. (fT§)) . the frequen- 
cies of lowest upper and lower polariton modes (m = 1 
is omitted in following discussion) are determined by 
n(f2]j t u) = w ex /f2L,u, and Eq. (J2TJ is rewritten as 



kmbc(^l,c/) 



1 + (^/u; ex ) 2 ' 



(25) 
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FIG. 2. (a) Frequencies of lower and upper polaritons (LP 
and UP) are plotted versus g/u cyL . Dashed curves are calcu- 
lated by Hamiltonian ©, while solid curves are obtained by 
Eq. (|19|) . (b) Dissipation rates of the polaritons are plotted. 
Dashed curves are calculated from Eq. (JTJ) with the RWA and 
ko/i^cx = 10 -2 . Solid curves are obtained from Eq. (|21[) (de- 
rived by the MBCs) by supposing A = 7.822 and I tuned for 
satisfying tan(oj cx i/c) — 1/A. Dotted line represents Ko/wex- 



If the light- matter coupling is not so strong Att(3 <C 1, we 
get Ql/u — w cx ± g, and kmbc(^7l,£/) is approximately 
a half of kq. Then, in the weak and normally strong 
coupling regimes, the resonance frequencies of cavity po- 
laritons and their dissipation rates certainly agree with 
the well-known ones. 

In Fig. rjl we plot (a) frequencies Ol,u of lower and 
upper polaritons (LP and UP) and (b) their dissipation 
rates as functions of gfc =W( . x / c /w x = y/4ftj3/2. For dashed 
curves, Qljj are calculated by diagonalizing Eq. at 
k = cj cx /c, and the dissipation rates are derived from 
Eq. ([1]) with K m =i(u) = k = 10~ 2 cj cx with the RWA 
on the system-reservoir coupling. For solid curves, fih,u 
and %bc(^l,[/) are determined by Eqs. ([lUf and (|2"T1) . 
respectively, for A = 7.822 corresponding to kq/lo ck — 
1CP 2 . As seen in Fig. [2] the two approaches give almost 
the same results in the weak and normally strong light- 
matter coupling regimes g <C w C x, where the frequencies 
are given as £2l,u — w cx ± g and the dissipation rates are 
Ko/2. In the ultrastrong light-matter coupling regime 
g > w cx , concerning the resonance frequencies the 
two approaches agree well with each other at least qual- 
itatively, while there is a discrepancy for the lower po- 
lariton frequency Ql (this is because Ql closes to the 
dissipation rate k ). However, as seen in Fig. [2{b), the 
dissipation rates obtained by the two approaches show 
qualitatively different behaviors. The solid curves can be 
well fitted by Eq. (|25|) [28| . while the dashed curves are 
expressed aa K f^ A = K,e/|% 0,111, 

where coefficient 



wl,u appearing in Eq. ([5]) is determined by the diagonal- 
ization [29(. In this way, in the ultrastrong light- matter 
coupling regime, Eq. ([TJ give incorrect dissipation rates 
compared with Eq. (f2~Tj) determined by the MBCs. 

In summary, based on the QED theory for dielectric 
media [HrlliJ , from the weak to ultrastrong light- matter 
coupling regime, the quantum Langevin equation (|24ap 
of discrete cavity polariton modes and the input-output 
relation ()24b[) are derived in the good cavity limit. They 
suggest that Eq. ((2]) is appropriate for the coupling be- 
tween the polariton system and the photonic reservoir in 
the good cavity case (30l |. This result is also consistent 
with Ref. ljl, 24 1, the discussion of the universe modes 
31 - 34 1 , and the Freshbach's projector approach 3|| 3|| , 
while empty cavities were mainly discussed. Although 
the counter-rotating terms such as a(ui)pj and p^a^ui) 
might appear in g eneral due to the overlap between the 
cavity modes [321 135.]. they are negligible thanks to the 
quality of the cavity, not to the RWA on system-reservoir 
coupling. As seen in Fig. HJb), the dissipation rates of 
the cavity polaritons derived from Eq. ([T]) with the RWA 
show the qualitatively different behavior compared to 
Eq. (J2TJ) derived by the QED theory (MBCs), whereas 
they agree well with each other in the weak and nor- 
mally strong light-matter coupling regimes. This can be 
understood in the sense of perturbation. If the light- 
matter coupling is stronger than the coupling with envi- 
ronment, we should diagonalize the cavity system first, 
then the coupling with environment is treated perturba- 
tively. Then, Eq. ((2|) i s ap propriate in general for the en- 
tire coupling regimes [37[ ■ If we want to discuss also the 
bad cavity case, we have to explicitly consider the MBCs 
(flTf . and the simplified expressions (fl) and ((2|) cannot 
be used. Whereas we simply supposed the cavity sys- 
tem depicted in Fig. [TJ specific cavity structures can also 
be considered based on the QED theory for dielectrics 
[ifl If we consider a loss or a gain inside the 

cavity by supposing the imaginary part of the dielectric 
function s p (uj), the Langevin equation (|24a[) would have 
another noise term expressed by Jn(0 < z < /,w). This 
work remains for the future. Further, whereas the sim- 
ple dielectric (polariton) medium was discussed in this 
letter, for cavity systems with an atom (nonlinear sys- 
tems) and for the Dicke model (the gauge invariance is 
broken) (38l - l4oj |. it is still open to dispute whether the 
system-reservoir coupling is generally expressed by the 
eigen-states of the cavity system as in Eq. ((2]). How- 
ever, the result in this letter suggests a universal policy: 
The system-reservoir coupling must be determined from a 
fundamental viewpoint after diagonalizing strongly cou- 
pled composite systems. Even for the superconducting 
circuits 11- 111, the system- reservoir coupling under the 
ultrastrong coupling between artificial atoms and res- 
onator modes should be determined by an appropriate 
microscopic description |4lj |. 
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Quantum fluctuation in homogeneous dielectric system 



In an homogeneous medium with dielectric function e(z, uS) = e(w), from the wave equation (|10[) . the vector potential 
is expressed as 



A + (z,u>) = A+(z,uj) + i+0,w), 



where At, and A< are forward and backward fields defined as 



A+(z,uj) = -fj, 



e ik(uj)(z-z') ^ 

dz' — tttt-, — ; — Jn(z',uj), 



— oo 

oo 



dz' 



i2k(uj) 

-\k{u)){z — 



i2k{u) 



-Jn(z',w). 



(26) 

(27a) 
(27b) 



Since the noise current density J-^(z,w) has the local correlation (commutablc for different positions z ^ z'), the 
forward field A > {z) and backward one A < (z') are commutablc for z < z' 



[A+(z,w),A<{z',w']\ =0 for z < z . 
The quantum fluctuation of the vector potential in the frequency domain is obtained as 

[A+(z, cA Az(z', = S(u - uj') MMe^M](«^)-i»W"*-/| 

AtteqcijjS \n{uj)\ 2 

When we define the spatial Fourier transform for wavenumber q > as 

1 



A+ - 



dz e- iqz A+(z,uj), 



the quantum fluctuation of this mode in the frequency domain is obtained as 

2Im[fc(w)] h Re[n(cj)] 



[A+( U ),Af(u/)] = 8(w- U ')8 qt f 

= S(UJ — u')5q t , 



q- k(uj)\ 2 AttEqCloS \n(uj)\' 



' q ' 47T£ ccj5 |n(w)| 2 
In the loss- less limit (Im[e(w)] — > 0), we get 



dz e i[q - k(u)]z + dz e l[q - k{ui) " ]z 



[i+M, Az(u')] = S(u - ui')8 q , q ,8{q - 



1 



q \-„--g-v~Ji - v~ -y~„ 2 £qCU jS n{u) ' 

Then, the equal-time quantum fluctuation of this mode is finally written as 



(28) 
(29) 

(30) 

(31a) 
(31b) 

(32) 



6q, q > 



duj 8(q — k(uj)) 



1 



1 



2£qooS n(u) 



2eocqS n{Q q ) ' 



(33a) 



where Q q = cq/n(fi q ). Therefore, in the loss- less dielectric media, the quantum fluctuation of the vector potential is 
modified by the factor of n(i? g )~ 1 . The fluctuations of the other electromagnetic fields are obtained as 



[£+ BZ] = 8„ 

[£>+,£)-,] 



hcq 



1 



M 2e Q Sn{n q f' 
hq 1 
2e cSn(Q q y 
hcq 



2e Q S 



n(ft q 



(33b) 
(33c) 
(33d) 
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The dependence of these fields on n(fl q ) is exactly the same as the one discussed in Ref. [17| . Then, the electromagnetic 
fields are sub-fluctuant or super-fluctuant in dielectrics compared to the case in vacuum. 
The forward and backward fields ([27)1 are rewritten as 



= J- » f , K Re[ ;y ] e i M^)]^ > (,,a;) 1 (34a) 

y 4neQCUjRe[n(uj)\b n{uj) 



A+ <M = \ A p r f , 1g Re[ ? ( T )] e- iRe[fcM] ^<(^^), (34b) 

< y 47T£ cwRe[n(w)]S n{ui) 



where operators a> and a< are defined as 



& < (z,u) = i. L "f°°f - e'"W»)]. / d*' e"<-*' ./„(*». (35b) 



/tu;Re[n(aj) 

They correspond to the annihilation operator of a photon in the dielectric medium, and the commutator is derived as 

[a > (z,cj),al{z , ,oj)} = [a K (z, u), (z', w)] = 8(u - w / )e- Im l*(«'0]|«-*'l ( 36 ) 

In the loss-less limit (Im[e(cj)] — > 0), they become position-independent and simply considered as the annihilation 
operator. 

Diagonalization of polariton system 

The Hamiltonian of polariton system is expressed in Eq. ([3]) . The vector potential and the electric field are expressed 
by a k as 



Ak = ^2^ks^ + ^> (37a) 



< hc\k\ 

The excitonic polarization and current density arc written in fc-space as 



A = H/™(a*-4)- ( 37b ) 



Pk = \l — s — ^ k k ^ ^ ' 



J fe = (-itd„W ' " LX (b fc - V_ k ) - 4npe uj cx 2 A k . (37d) 



The equations of motion are derived as 

d_ 

di 
d_ 

di 

The equations of motion of the macroscopic fields are obtained as 



a k = -ic\k\a k + 9k(h - b- k ) - i2D k (a k + a'_ k ), (38a) 
bk = -iu eK b k - gk(ak + al k ). (38b) 



^i fc = -E k , (39a) 
at 

^-E k = c 2 k 2 A k - -J fe , (39b) 
at £o 

4 A = Jk, (39c) 
at 

— J k = -uj cx 2 P k - 47r/3e w ox 2 1 T:ifc- (39d) 
at at 
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The first two equations correspond to the Maxwell equations. Then, the wave equations are obtained as 

1 d 2 - 
^ Ak + ^2 ~Qt2 Ak = ( "° Jfe ' 



, 2 ; I d 2 - 



'dt 2 



The equation of motion of the polarization is rewritten as 

d 2 



dt 2 



P k = -W ex Pk + 47T/3iToWex E k . 



Then, by Fourier transforming to the frequency domain, we get the dielectric function as 



21,2 



c 2 k 



1 + 



4nf3uj c 



2 , ,2 



W w cx - (cj + i0+) 2 cj cx - (w + i0+) 2 ' 

Here, ix> ox l is the frequency of longitudinal excitation satisfying 

2 

Att/3 = - 1. 



(40a) 
(40b) 

(41) 

(42) 

(43) 



The Hamiltonian can be diagonalized by the Bogoliubov transformation. The annihilation operator of lower (j = L) 
or upper (j — U) polariton is represented as 



Pj,k = Wjk&k + Xj k bk + VjkoLk + z jk V_ k . 

Then, the coefficients are determined by the following eigen-value problem: 

(c\k\ + 2D k -ig k 
igk w cx 
'-Ih, -ig k -c|, 
\ -igk 

The four eigen-states corresponds to pj_ k and k , whose eigen-frequencies are respectively Wj, k and —u)j t k : 



-2D k 




( W 3,k\ 




( w oA 


—ig k 





3->k 






k\ - 2D k 




Vj-k 


y 3 ,k 


igk 




\ z 3,k) 




\ z j,kJ 



(44) 



(45) 



u L /u,k = \ 1 + 4tt^ - 



c 2 |fc| : 



1 + 4tt/3 + 



c 2 |fc| 2 



The eigen vectors are derived as 



( 

XL,k 
UL,k 
\ZL,k ) 



w cx - 



47T/3 



-1/2 



4c 2 \k\ 



Lj L +c\k\ 



1/2' 



1/2 



(46) 



W ox 2 J 20Jex V 

l -^(i-^) ; 



(47) 







Xu.k 




UU,k 




\ z U,k J 





1 - 



UJ C 



+ 47T/3 



-1/2 



V 



iV^( 

— bJu 2 . 

"ex 2 



wt/+c|fc| 

2u>ex 



LJ[/ — c|fc| 

2w cx 



(48) 



We have two eigen-frequencies uij. k for a given wavenumber k. Inversely, we get one allowed wavenumber k(u>) for a 
given frequency w. The relation of them are expressed by the dielectric function as 



c 2 k 2 
^ 3 ,k 2 
c 2 k(uj) 2 

OJ 2 



e(co). 



(49a) 
(49b) 



The Hamiltonian is rewritten as 



E E ^3,w],kP3,k 

j=L,U k— — oo 
poo p -, 

= J du ha; \pl (w)p> (w) + p], (uj)p< (u) , 
where the forward and backward polariton operators (p> and p < ) are defined in the frequency domain as 



L 



2irv g (uj) 



[0(uj - UexL,)Pu,±k(_u>) + 0{u ~ Wex)Pi,±fc(w)] 



(50a) 
(50b) 

(51) 



We choose the phase difference between lower and upper polaritons by the phase of the eigen- vectors shown above for 
deriving a simple expression below. The group velocity is expressed as 



dk 



- 1 



- 1 + 47T/3 



By using the polariton operator, the original photon and excitation operators are expressed as 

Ofc = E ( w j,kPj,k-yj,kPj,- k )i 

j=L,U 

h k = E ( X lkPj,k- Z jlk pt_ k ). 
j=L,U 

Then, the vector potential at position z is represented as 

1 oo 

Mz) = j= E eikzA * 



fc— — c 



E 

k— — oo 
oo 

E 

k— — oc 



2e c\k\SL 



(a k + a_ k )e 



E _ y*j-k)Pj,k + {w 3 ,-k - yj,k)p],- k 



2e c\k\SL.^ u < 



'ikz 



The positive-frequency and forward component is then written as 

oo / ^ 

^(*) = E ' 



k=0 V " 1 1 j=L,U 

By rewriting the polariton operator in the frequency domain, we get 



4tw= E E 

j=L,U fe=0 



E E 

j=L,U k— — oc 
i / 



h 


2 


ck 


2e ckSL 

V 






2 

+ 47T/3 



PjM 



I hv g {uj 3k ) - kz 

2e c 2 kSL Pjke 



n 



47r£ cwn(a;)5' 

This has exactly the same form as Eq. (1541 in the loss-less limit. 



P>(")e ife(u 



(52) 

(53a) 
(53b) 



(54a) 
(54b) 
(54c) 



(55) 



(56a) 

(56b) 
(56c) 
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Solution in cavity system by Green's function approach 



Let's derive the Green's function G(z,z',u) satisfying 

where the dielectric function is expressed as 

e(z,u) = r](w)d(z) 



G(z,z',ui) = S(z - z'). 



1 z < 

e p (w) < z < I 



(57) 



(58) 



Obeying the recipe in Ref. [42|, section 2.4], the Green's function Gy(z, z' ,lu) (z in region i and z' in region j) is 
obtained as 



Gu(z,z',oj) 
G 2 i(z,z',uj) 
Gl2(z, z',(J) 
G 22 (z,z',uj) 



1 f e i( M /c)|g-»'l +e -^h)*r 2 x{u)e- i{u/c)z '\ 
\2{(jjIc) L J 

sm[k p (u;)(l ~ z)}^^/^', 



i2(w/c) 
1 

i2(w/c) 
1 

i2/c p (w) 



(59a) 
(59b) 
(59c) 



1 <r e ifcp 

i2/c p (w) I 



; p (w)|*-«'| _ -ik p (uj)(z-l) -ik p (uj)(z' -I) 



sm[kp(uj)(l - z)\ sin[fc p (cj)(2 - z')] \ (59d) 



[1 - LA(w)] 8in[jfcp(w)Z] + m p (w) cos[fc p (w)Z] pV A p 

e ifcp(o>)|*-*'| _ c ifc„Mz I 1 Z iyl M ~ "pHI sin[fc p (q;)(^ - z')] 



i2fc p (o;) [1 — LA(w)] sin[fep(cj)Z] + inp(cj) cos[fc p (w)Z] 

-lfc p ( 2 -0 [1 - sin[fc p (q;)z / ] + m p (u) cos[fc p (q;)z / ] 

[1 — sin[fc p (w)Z] + in p (cj) cos[fc p (w)Z] 



?2lM = 
*2l(<*>) = 



[1 + L/l(w)] sin[/c p (w)/] — m p (w) cos[fc p (w)Z] 
[1 — L4(w)] sin[fe p (o;)2] + irt p (w) cos[fc p (o;)Z] 
2 

[1 — sin[/c p (o;)/] + in p (u;) cos[fc p (o;)Z] 

2n. p (w) 



[1 — L4(w)] sin[fc p (o;)£] + in p (uj) cos[fc p (o;)Z] 
The derivation is as follows. When a source exists in region 1, the Green's function can be supposed as follows: 

Uji{Z,Z ,UJ — — — — , 

i2(w/c) 

g 11 (z,z',uj) = e i ^^ z - z ^+e- i ^ c >B 11 (z'), 
Q 21 {z,z\u:)=^ k ^ z -^F 21 {z')+e- ik ^ z -^B 21 {z'). 

We get a boundary condition at z = I: 

F 21 {z') + B 21 (z') = 0. 

At z = 0, we also get 

e -i(u/c)z' + Bn ( z ') = e - ife p(") z F 2 i(z') + e ik " {uj)l B 21 (z'), 
e -i(«/o)«' _ _ np ( w ) [e- ifc p^ z F 21 (z') - e ifc ^B 2 i(z')l = -LA(w) [e" 1 ^/ ) 2 ' + Bn^) 



(59e) 

(60a) 
(60b) 
(60c) 

(61) 

(62a) 
(62b) 

(63a) 

(63b) 
(63c) 
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The third condition is obtained by the boundary condition (jTTbJ) or integrating Eq. (|61l) . Solving them, we get Gn 
and G21. 

When a source exists at region 2, we can suppose 

Gi2{z,z,ui) = - . 2fcpM , (64) 

G 12 (z,z',u) = e-^/^B^z'), (65a) 
g 22 (z,z , ,uj) = e ifc p(")l z - z 'l + e^"^ 2 F 22 (z') + e- ik ^ )(z - l) B 22 (z'). (65b) 



At z = I, we get 



Further, at z = 



i*p(«)(i-*') + e^^'-FaaOO + B 22 (z') = 0. (66a) 



e 



M")*' + F22 ( z ') + jW l B 22 {z') = B 12 (z'), (66b) 



S 12 (z') + np(«) [-e"^ 2 ' + F 22 (z') - e ife ^B 22 (z')] = LA(w)Bi2(z')- (66c) 

Then, G 22 and Gi 2 are obtained. 

By using the Green's function, the vector potential is obtained in region 1 as 



Af(z, u) = jWJ'A^ (z, u) + e-^' c >hiAi> (0, w) 





dz' e-^^X'-^JN^') +Mo / dz' Gi 2 (z, z', u) J N (z'). (67) 



i2(w/c) J z ' ' Jq 

where we define the incoming field as 



A+ >( ^) = - f dz' e- Wj N (z'). 

i2(w/c) ./_«, 



(68) 



Here, as discussed in Ref. 22J, we can focus only on the incoming and outgoing fields at the boundary z = 0, because 
the contribution such as the third term in Eq. (|67[) can be neglected in the loss-less limit (if there is no dissipation, 
we need not consider the noise operator for the propagation in vacuum). Further, the last term can also be neglected 
in the loss-less limit, then the vector potential in region 1 can be expressed as 

A+{z,uj) = e^W^^j+e-'W^i^u). (69) 

The first term corresponds to the incoming field A^^oj) = ^^(O, cj), and the second term is the outgoing field 
represented as 

A+ < (u)=r 2 iA+ > (0,u). (70) 
On the other hand, the vector potential in region 2 is obtained in the loss-less limit as 

i+(z,u;) =t 21 Af > (u;)sin[k p (uj)(l-z)]. (71) 

Then, we get 

A+{u:)=t 21 A+ > {uj). (72) 

These fields A^iin), Aj^w), A^ioj) certainly satisfy the Maxwell's boundary conditions. 
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